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Abstract 

In this paper we consider the inverse problem of determining a rigid 
inclusion inside a thin plate by applying a couple field at the boundary 
and by measuring the induced transversal displacement and its normal 
derivative at the boundary of the plate. The plate is made by non- 
homogeneous linearly elastic material belonging to a general class of 
anisotropy. For this severely ill-posed problem, under suitable a priori 
regularity assumptions on the boundary of the inclusion, we prove a 
stability estimate of log-log type. 
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1 Introduction 

The present study continues a line of research on Non-Destructive Tech- 
niques (NDT) in the theory of elastic plates. Specifically, the paper deals 
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with the determination of a rigid inclusion embedded in elastic plates by 
measurements taken at the boundary. This inverse problem arises in damage 
assessment of plates which are possibly defective due to the presence of in- 
terior rigid inclusions induced during the manufacturing process. The most 
simple and common NDT is the visual inspection of the specimen. However, 
the visual inspection is unable to detect a rigid inclusion which is embedded 
in a plate or which is indistinguishable from the surrounding material. Other 
conventional NDTs, such as those based on thermal or ultrasonic analysis, X- 
rays methods and others, are local by nature. To be effective, these methods 
require the vicinity of the inclusion be known a priori and readily accessible 
for testing, whereas in many practical situations the boundary of the plate is 
the only accessible region. Vibration response-based identification methods 
were also recently developed to overcome these difficulties, see, for example. 



The basic idea of most NDTs is that the inclusion changes the response of 
a plate and that the defect can be determined by comparing the response of 
a possibly defective specimen with a reference specimen, i.e. a plate without 
inclusion. In the case of inclusion made by (different) elastic material, in 
|M-R,o-Ve3j . |M-R,o-Ve5] constructive upper and lower estimates of the area of 
the inclusion in terms of the difference between the work exerted in deforming 
the defective and the reference plate by applying a given couple field at the 
boundary were determined. These results give no indication on the position 
of the domain occupied by the inclusion and, unfortunately, they still have 
not been extended to the case of a rigid inclusion. In this paper we tackle 
these two issues and we establish a stability result for the inverse problem of 
determining a rigid inclusion in an anisotropic elastic plate by one boundary 
measurement. 

From the mathematical point of view, [Fi], |Guj . we shall work within the 
Kirchhoff-Love theory of thin, elastic anisotropic plates under infinitesimal 
deformations. Let Q denote the middle plane of the plate and let h be its 
constant thickness. We assume that is a bounded domain of of class 
C^'^. The rigid inclusion D is modelled as an open simply connected domain 
compactly contained in Q. The transversal displacement w G H'^{Q) of the 
plate satisfies the following mixed boundary value problem 



El- 




mQ\D, 



on dQ 



(1.3) 



w)n-T),s= (Mr) 



on dQ 





in D 



dw" Qui' 

dn dn ' 
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coupled with the equilibrium conditions for the rigid inclusion D 

(1.6) / {dw{FV^w)-n + {{FV^w)n-T),s)g-{{fV^w)n-n)g^n = 0, 

JdD 

for every g E A, 

where A denotes the space of affine functions. In the above equations, n and 
r are the unit outer normal and the unit tangent vector toQ\D, respectively, 
and we have defined w'^ = w\q\jj and = wljy. Moreover, Mr, Mn are the 

twisting and bending components of the assigned couple field M, respectively. 
The plate tensor P is given by P = f^C, where C is the elasticity tensor 
describing the response of the material. 

Given any M G H~^{dfl,M.'^), satisfying the compatibility conditions 
Jg^Mi = 0, for i = 1,2, and for P bounded and strongly convex, problem 
(ll.ll) - fll.6p admits a solution w G H^{Q), which is uniquely determined up 
to addition of an affine function. 

The uniqueness issue for the inverse problem under consideration is the 
following. 

Given an open portion S of dQ and given two solutions Wi to (11. ip - (II. 6p 
for D = Di, i = 1,2, satisfying 

(1.7) Wi = W2, on S, 



dwi dwo „ 
1.8 —1 = —^^ on E, 

on on 

does Di = D2 hold? 

Let us notice that, by (II. 2p . (II. 3p . (II. 7p . (II. 8p . wi and W2 assume the same 
Cauchy data on S. Under the a priori assumption of C^'^ regularity of the 
boundary of the inclusion, a positive answer to the uniqueness issue has been 
given in |M-Ro2j . see also |M-Ro-Ve2 ] for a result of unique determination of 
a cavity in an elastic plate by two boundary measurements. 

In order to face the stability issue, following the approach used for the 
case of rigid inclusions in elastic bodies in |M-Rol] . we found it convenient 
to replace each solution Wi introduced above with Vi = Wi — gi, where gi 
is the affine function which coincides with Wi on dDi, i = 1,2. By this 
approach, maintaining the same letter to denote the solution, we rephrase 
the equilibrium problem (ll.ip - (ll.5p in terms of the following mixed boundary 
value problem with homogeneous Dirichlet conditions on the boundary of the 
rigid inclusion 



3 



(1.9) 

(1.10) 

(1.11) 

(1.12) 

(1.13) 



' div(div(PV2u;)) = 0, 
{FV'^w)n-n = -M„, 

< dw(FV^w)-n+{{FV'^w)n-T)„={Mr) 
w = 0, 



mn\D, 
on dQ, 
on dQ, 
on dD, 
on dD, 



dw g 

^ dn 



which has a unique solution w G H^{Q \ D). 

On the other hand, it is clear that the arbitrariness of this normaliza- 
tion, related to the fact that gi is unknown, i = 1,2, leads to the following 
formulation of the stability issue. 

Given two solutions Wi to f ll.9p -( lLT3|) . (11. 6p when D = Di, i = 1,2, 
satisfying, for some e > 0, 



to evaluate the rate at which the Hausdorff distance between Di and D2 tends 
to zero as e tends to zero. 

In the present paper, assuming C^'^ regularity of dD, we prove the fol- 
lowing constructive stability estimate of log-log type 



where C, r], C > and < 77 < 1, are constants only depending on the a 
priori data, see Theorem 13.11 for a precise statement. 

The methods used to prove (11.151) are based essentially on quantitative 
estimates of unique continuation from Cauchy data for a solution of the mixed 
problem f ll.9p -( TTTT^ . (II. 6p . These estimates involve a cornerstone result of 
unique continuation, namely the Three Spheres Inequality (14. 2 p for solutions 
to the plate equation (II. 9p . which has been determined in [M-Ro-Ve5] under 
the very general assumption that the elastic material of the plate obeys the 
dichotomy condition (I2.30ap -f l2.30bp . 

The logarithmic character of the stability estimate (ll.lSp seems difficult to 
improve. First, the Cauchy problem up to the boundary is severely ill-posed 
and, even in the simpler context of the electrical impedance tomography 
which involves a second order elliptic equation instead of a fourth order one, 
the corresponding stability estimate with a single logarithm is the best possi- 
ble result, see |Al-B-Ro-Vej and also |Al-R-Ro-Vej for a general discussion on 
the ill-posedness of the Cauchy problem. We also quote |DiC-Rj for examples 
of exponential instability for the inverse inclusion problem in a conducting 
body. In addition, a further difficulty arises in our analysis. It is due to 



(1.14) 




(1.15) 



(i^(A,D2) <C(log|loge|) 
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the lack of quantitative estimates of the strong unique continuation property 
at the boundary in the form of either Three Spheres Inequality or Doubling 
Inequality. It has been shown in |Al-B-Ro-Ve] that this is a key ingredient 
in proving that the stabihty estimate for the corresponding inverse prob- 
lem with unknown boundaries in the conductivity context is not worse than 
logarithm. This mathematical tool is not currently available for the plate op- 
erator, even in the simplest case of isotropic material, and this is the reason 
for the presence of a double logarithm in the stability estimate f ll.lSp . We re- 
fer to |M-Rol] for a discussion on the analogous problem in the determination 
of a rigid inclusion in an isotropic elastic body by boundary measurements. 
Finally, as remarked in |M-Ro-Ve4] , it seems hopeless the possibility that so- 
lutions to f ll.ip can satisfy even a strong unique continuation property at the 
interior, without any a priori assumption on the anisotropy of the material, 
see also [Alij. Regarding this point, our dichotomy condition fl2.30ap -( l2.30bp 
basically contains the same assumptions under which the unique continuation 
property holds for a fourth order elliptic equation in two variables. 

The paper is organized as follows. Some notation is presented in Section 
2. In Section 3 we state two auxiliary propositions concerning the estimate 
of continuation from the interior (Proposition 13. 2p and from Cauchy data 
(Proposition 13. 3p . and we give the proof of the main Theorem 13.11 Section 
4 contains the proofs of Proposition 13.21 and Proposition 13.31 Proofs of regu- 
larity estimates of the solution to the mixed problem fll.9p - fll.13p . (11. 6p are 
presented in Section 5 and, in part, in the Appendix. 

2 Notation 

Let P = (xi(P),X2(P)) be a point of M^. We shall denote by Br{P) the 
disk in of radius r and center P and by Ra^b{P) the rectangle of center 
P and sides parallel to the coordinate axes, of length 2a and 26, namely 
Ra,b{P) = {x = {xi,X2) I l^i — xi{P)\ < a, \x2 — X2{P)\ < b} . To simplify 
the notation, we shall denote Br = Br{0), Rafi = Ra,b{0). 
Given a bounded domain Q in we shall denote 

(2.1) Qp = {x eQ\ dist(x, dQ) > p}. 

When representing locally a boundary as a graph, we use the following defi- 
nition. 

Definition 2.1. (C'^'^ regularity) Let f2 be a bounded domain in M^. Given 
A; G N, we say that a portion S of dfl is of class C^'^ with constants po, 
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Mq > 0, if, for any P E S, there exists a rigid transformation of coordinates 
under which we have P = and 



AIq ''^^ AIq 

where -0 is a C'''^ function on (^—-^, satisfying 

^(0) = 0, 
^/;'(0) = 0, when k > 1, 
IIV'ILfc,ir-«L «L^ < ^Vo- 

\ Mo' Mo J 

When = we also say that S is of Lipschitz class with constants Pq, Mq. 

Remark 2.2. We use the convention to normahze all norms in such a way that 
their terms are dimensionally homogeneous with the L°° norm and coincide 
with the standard definition when the dimensional parameter equals one. For 
instance, the norm appearing above is meant as follows 

fc+i 

V «0'«0y ^ V Mo'MoJ 

where i/j^^^ is the z-th derivative with respect to the xi variable. Similarly, 
denoting by V^u the vector which components are the derivatives of order i 
of a function u defined in Q, we denote 

fc+i 

ikiic"=.i(f^) = y^po'iiv''»iu°°(Q), 

i=0 



.1=0 

and so on for boundary and trace norms such as 

Notice also that, when Q = B^, then Q satisfies Definition 12.11 with 
Po = R, Mq = 2 and therefore, for instance, we have 



\i=0 
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Given a bounded domain Q in such that dfl is of class C'^'^, with 
k > 1, we consider as positive the orientation of the boundary induced by 
the outer unit normal n in the following sense. Given a point P G dQ, let 
us denote by r = t{P) the unit tangent at the boundary in P obtained by 
applying to n a counterclockwise rotation of angle |, that is 

(2.2) T = 63 X n, 

where x denotes the vector product in M^, {ei, 62} is the canonical basis in 
and — ei x 62- 

Given any connected component C of dQ and fixed a point P & C, 
let us define as positive the orientation of C associated to an arclength 
parametrization ip{s) = {xi{s),X2{s)), s e [0,/(C)], such that (/?(0) = P 
and (f'{s) = t{(p{s)), where 1{C) denotes the length of C. 

Throughout the paper, we denote by diU, dsU, and dnU the derivatives of 
a function u with respect to the Xi variable, to the arclength s and to the 
normal direction n, respectively, and similarly for higher order derivatives. 

We denote by the space of 2 x 2 real valued matrices and by C{X, Y) 
the space of bounded linear operators between Banach spaces X and Y. 

For every 2x2 matrices A, B and for every L e >C(M^,M^), we use the 
following notation: 



(2.3) (I-'^)ii — LijkiAki, 

(2.4) A-B = AijBij, 

(2.5) \A\^{A-A)^, 

(2.6) A'y"^ = ^ (A + A*) , 



where A* denotes the transpose of the matrix A. Notice that here and in the 
sequel summation over repeated indexes is implied. 

Finally, let us introduce the linear space of the affine functions on R^ 

A — {g{xi, X2) — axi + bx2 + c, a,b,c& R}. 

2.1 A priori information 

i) A priori information on the domain. 
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Let us consider a thin plate Q x [— | , |] with middle surface represented by 
a bounded domain Q in and having uniform thickness h, h « diam(Q). 
We shall assume that, given po, -^i > 0, 

(2.7) \n\ < MipI 

where \Q\ denotes the Lebesgue measure of Q. We shall also assume that fl 
contains an open simply connected rigid inclusion D such that 

(2.8) dist{D,dn)> po. 

Moreover, we denote by S an open portion within dfl representing the part 
of the boundary where measurements are taken. 

Concerning the regularity of the boundaries, given Mq > 0, we assume 
that 



(2.9) dQ is of class C ' with constants po, M( 



0; 



(2.10) E is of class C ' with constants po, Mq. 

(2.11) dD is of class C^'^ with constants po, Mq. 

Moreover, we shall assume that for some Pq G E and some 5o, < 5o < 1, 

(2.12) aoni?^,,„(Po)cS, 

and that 

(2.13) \E\ < {1 - 6o)\dn\. 

a) Assumptions about the boundary data. 
On the Neumann data M we assume that 

(2.14) M e L^{dn,R''), (M„,(M,)„) ^0, 

(2.15) supp(M) CC E, 
the (obvious) compatibility condition 

(2.16) / Mi^O, i = l,2, 

Jan 



and that, for a given constant F > 0, 

\\M\\L^dn,R^) ^ ^ 

Hi) Assumptions about the elasticity tensor. 

Let us assume that the plate is made of nonhomogeneous hnear elastic 
material with plate tensor 

(2.18) ^=1^^' 

where the elasticity tensor C{x) G C(NP, M^) has cartesian components Cijki 
which satisfy the following symmetry conditions 

(2.19) CijM = Ckiij = Ckiji k,l = 1, 2, a.e. in Q. 
We recall that fl2.19p are equivalent to 

(2.20) CA = CA'y"^, 



(2.21) 
(2.22) 



CA is symmetric, 
CA - B = CB ■ A, 



for every 2x2 matrices A, B. 

In order to simplify the presentation, we shall assume that the tensor C 
is defined in all of M^. 

Condition f l2.19p implies that instead of 16 coefficients we actually deal 
with 6 coefficients and we denote 



(2.23) 

and 

(2.24) 



Ciiii 


= ^0, 


C1122 = 


C2211 = Bq 




Clll2 


= C'1121 


= C'1211 


= C'2111 = 


Co 


< C2212 


= ^72221 


= (^1222 


= C'2122 = 


Do 


C1212 


= ^71221 


= (^21 12 


= C'2121 = 




C2222 


= F,, 









ao = Ao, ai = 4Co, aa = 2So + AEq, = ADq, 04 = Fq. 
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Let S{x) be the following 7x7 matrix 



(2.25) 



S{x) 







ai 


a2 


"3 


04 





\ 







ao 


ai 


^2 


as 


04 













ao 


ai 


a2 


as 


a4 




4ao 


3ai 


2a2 


as 
















4ao 


3ai 


2a2 


as 
















4ao 


3ai 


2a2 


«3 





V 











4ao 


3ai 


2a2 


as J 



and 
(2.26) 



V{x) = —\detS{x)\. 
ao 



On the elasticity tensor C we make the following assumptions: 
I) Regularity 

(2.27) C e C^•^(M^£(M^M2)), 

with 



(2.28) 



2 2 

i,j,k,l=l m=0 



where M is a positive constant; 

II) Ellipticity (strong convexity) There exists 7 > such that 



(2.29) CA-v4>7|A|2, inM^ 

for every 2x2 symmetric matrix A. 
Ill) Dichotomy condition 

(2.30a) either T){x) > 0, for every x E M^, 

(2.30b) or V{x) = 0, for every x E ]R^ 

where V{x) is defined by f l2.26p . 

Remark 2.3. Whenever fl2.30ap holds we denote 

(2.31) Si=mmV. 

We emphasize that, in all the following statements, whenever a constant is 
said to depend on 61 (among other quantities) it is understood that such 
dependence occurs only when fl2.30ap holds. 
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We shall refer to the set of constants Mq, Mi, Sq, F, 7, M, 5i as the a 
priori data. The dependence on the thickness parameter h will be omitted. 

In the sequel we shall consider the following boundary value problem of 
mixed type 



(2.32) 
(2.33) 
(2.34) 
(2.35) 
(2.36) 



f div(div(PV2«;)) = 0, 
(FV'^w)n-n = -M„, 

dw{FV^w) ■ n + {{FV'^w)n ■ r)„ = (M^)„ , 
w = 0, 

dw ri 

dn ~ ^' 



mn\D, 
on dQ, 
on dQ, 
on dD, 
on dD, 



coupled with the equilibrium conditions for the rigid inclusion D 
(2.37) / (div(PV^'«;) ■ n + ((PV^ w)n ■ t),s 

JdD 

for every g E A. 



By standard variational arguments, it is easy to see that problem fl2.32p - 
fl2.37p admits a unique solution w G H^{Q \ D) such that 

(2.38) < Cpl\\M\\^-^,^,^^^.y 

where C > only depends on 7, Mq, and Mi. 



3 Statement and proof of the main result 

Here and in the sequel we shall denote by G the connected component of 
VL \ {Di U D2) such that S C dG. 

Theorem 3.1 (Stability result). Let Q be a bounded domain in satisfying 
(12. 7p and (12. 9p . Let Di, i = 1,2, be two simply connected open subsets of Q 
satisfying (12. 8 p and (12.111) . Moreover, let S be an open portion of dQ satisfy- 
ing f l2J0|) . (EUD and f l27[3|) . Let M e L'^{dn,R'^) satisfy fimj) - fl2Tn) and 
/et i/ie plate tensor F given by (I2.18P satisfy the symmetry conditions (I2.19p . 
the regularity condition (12.281) . the strong convexity condition (I2.29P and the 
dichotomy condition. Let wi G H^{Q \ D^) be the solution to (I2.32p - (l2.37p . 
when D = Di, i = 1,2. If, given e > 0, we have 



(3.1) 



mm 



\wi -W2- 5^1^2(2) + PO 



d_ 

dn 



>i -W2- g) 



< e, 
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then we have 

(3.2) dn{dD,,dD2)<poUJ 
and 

(3.3) dn(D;,D~2) < pou ^ 

where u is an increasing continuous function on [0, oo) which satisfies 

(3.4) uj{t) < C(log I logtl)"'', for every t, < t < e"\ 

and C,rj,C>0,0<r]<l, are constants only depending on the a priori 
data. 

The proof of Theorem 13.11 is obtained from the following sequence of 
Propositions. 

Proposition 3.2 (Lipschitz Propagation of Smallness). Let Q be a bounded 
domain in satisfying (12. 7p and (12. 9p . Let D be an open simply connected 
subset of Q satisfying (12.81) . (12. lip . Let w G H'^{Q \ D) be the solution 
to (I2:32|) - (12371 ) ■ where the plate tensor P given by fl238D satisfies (ICT]) . 
(I2.28p . (I2.29P and the dichotomy condition. Let the couple field M satisfy 

dZHD-dlllD. 

There exists s > 1, only depending on 7, M, 5i, Mq and 5q, such that for 
every p > and every x E {fl \ D)sp, we have 



(3.5) / 

JB 









exp 





\V^w?> ^ ^llMf 1 



Bp{x) 

where A > 0, B > and C > only depend on 7, M, 61, Mq, Mi, 6q and 
F. 

Proposition 3.3 (Stability Estimate of Continuation from Cauchy Data). 
Let the hypotheses of Theorem \3.1\ be satisfied. We have 

(3.6) / |VVr < Poll^^f -1 ^ 



D2\Dl " \^pg||M||^„i 
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(3.7) 



\V'w,\' < pI\\m\\ 



Di\D2 



CO 



pUm\ 



where u is an increasing continuous function on [0, oo) which satisfies 

(3.8) uj{t) < C{\og I logt|)"^, for every t < e^^, 

with C > only depending on 7, M, 61, Mq, Mi and 60. 

If we assume, in addition, that there exist L > and Pq, < po < Po, 
such that dG is of Lipschitz class with constants po, L, then fl3.6p -f l3l7|) hold 
with u given by 

(3.9) u{t) < C\ logtl'", for every t < 1, 

where a > and C > only depend on 7, M, 61, Mq, Mi, 6q, L and ^. 

Proof of Theorem 15*. i[ Let us denote, for simplicity, d = d^^i^dDi, dD-^). Let 
us see that, if > is such that 



(3.10) 

then we have 
(3.11) 



IV72 |2 ^ ^ 
_|V Wi\ < — , 

D2\Di Po 



_|V U)2\ < — , 
Di\D2 Po 



d < Cpo 



log 



'Cpl\\M\ 



where B > and C > only depend on 7, M, ^i, Mq, Mi, Sq and F. 

We may assume, with no loss of generality, that there exists xo G dDi 
such that dist(xo,5D2) = d. Let us distinguish two cases: 

i) Bdixo) C D2] 

ii) Bdixo) n = 0. 

In case i), by the regularity assumptions made on dDi, there exists Xi G 
D2 \ Di such that Btd{xi) C 1^2 \ -Di, with t 

By (I3.10p and by Proposition 13.21 with p = 



— , we have 



(3.12) 



Cpt 



exp 



m\ 



where A > 0, B > and C > only depend on 7, M, 5i, Mq, Mi, Sq and F. 
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By (KWf we easily find fIXTT]) . 

Case ii) can be treated similarly by substituting wi with W2. 
Hence, by Proposition 13.31 and assuming e < Po\\M\\ 1 we ob- 

tain 

(3.13) d<Cpo|log 

where B > and C > only depend on 7, M, ^i, Mq, Mi, Sq and F. 

Thus we have obtained a stability estimate of log-log-log type. 

In order to prove an analogous estimate for the Hausdorff distance be- 
tween Di and D2, let us now set d = d-^i^Di, D2) and let us assume, with no 
loss of generality, that there exists Xq G Di such that dist(xo,-D2) = d. If 
Bd^xo) C Di, then i?d(xo) C Di \ -D2 and fl3.12p follows with t replaced by 1. 
If, otherwise, Bd^xo) ^ Di, then dist(xo, i9i?i) < rf, and we can distinguish 
two cases: 

i) dist(a;o,aDi) > f, 

ii) dist(a;o,ai^i) < f. 

When i) holds, then Bd{xQ) G Di \ D2 and again (13.121) follows with t 

replaced by |. When ii) holds, there exists yo G dDi such that \yo — xo\ < |. 
Therefore there exists yi G Di such that Btd{yi) G Di\ D2, and (I3.12p 

follows with t replaced by |. From (I3.12p . arguing as above, we obtain (I3.13P 
for d = dn(D[,D2). 

Next, by this rough estimate, we can apply Proposition 3.6 in |Al-B-Ro-Vej 
which ensure that we can find eg > 0, only depending on 7, M, Mq, Mi, 
^0 and F, such that if e < eo then dG is of Lipschitz class with constants po, 
L, with L and ^ only depending on Mq. By the second part of Proposition 
K?X the log- log type estimates (O, (Q, (El follow. □ 



log 



log- 



pI\\m\\^^ 



5(9n,IR2) 



4 Proofs of Propositions 13.21 and 13.3 



We need to premise some auxiliary results. The first one is the basic tool of 
our approach, the three spheres inequality. 

Proposition 4.1 (Three Spheres Inequality ( |M-Ro-Ve5] . Proposition 5.1)). 

Let Q be a domain in M^, and let the plate tensor P given by (I2.18P satisfy 
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fl2.19p . fl2.28p . f l2.29p and the dichotomy condition. Let w G H'^iVt) he a weak 
solution to the equation 

(4.1) div(div(PV^«;)) = 0, m VL. 

For every ri, r2, r3, r, < ri < r2 < r^ <r, and for every x & fir we have 

(4.2) / \v'w\'<c( [ l^M'] ( I \^M^ , 

JBr^ix) \JBr-,{x) J \JBr.^{x) J 

where C > and 6, < 6 < 1, only depend on 7, M, 61, ^ and ^. 
Let us define, for p < po, 

(4.3) = {x G ^] I dist(x, an) < p}, 
tiiat isUf = n\ 

The following three lemmas state regularity estimates. 

Lemma 4.2 (Global regularity for the mixed problem). Let Q be a 

bounded domain in satisfying (12. 7p and (12.91) . Let D be a simply con- 
nected open subset of Vt satisfying fl2.8p and 02.111) . Let M e H^i^dVl,^^) 
satisfy (I2.16p . Let the plate tensor P be defined by (12.180 and satisfying 
fl219|) . ([228]), (12:291) . Let w E H^{n\D) be the solution to (12321) - fOTj) . 

We have 

where C > only depends on Mq, Mi, 7, M. 

Lemma 4.3 {H'^ regularity up to the boundary of the rigid inclusion). Let 
Q be a bounded domain in satisfying (12.71) and (12. 9 p . Let D be a simply 
connected open subs et ofVL satisfying and fl2lTD Let M e H-2{dn,R^) 
satisfy (I2.16p . Let the plate tensor P be defined by (12.180 and satisfying 
i^M, (^M- Let we H\n\D) be the solution to fl^ - fl07p . 

We have 

where C > only depends on Mq, Mi, 7, M. 
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Lemma 4.4 (C^'" regularity up to the boundary of the rigid inclusion). 

Let Q be a bounded domain in satisfying (12. 7p and (12. 9p . Let D be a 
simply connected open subset of Q satisfying (12. 8 p and (12. lip . Let M G 
H~2(^dQ,M.'^) satisfy (I2.16P and let the plate tensor F given by (12.181) satisfy 
fl2J9|) . (12:281) and ([229]). Lei w e H^{n\D) be the solution to fl2:32|) - flOTp . 

We have 

(4.6) lliyll ^ ,— > PON < Cpnll^ll 1. for every a < 1, 

^ ' " "c2.«((n\D)\w-r) — '^011 ii_f/-^(9n,iR2)' ^ ' 

where C > on/?/ depends on Mq, Mi, j, M , a. 

Proof of Lemma \4-4\ The above Lemma follows by Lemma and by Sobolev 



embedding theorem, [Ad] . □ 

Remark 4.5. We have stated the above regularity lemmas under the regularity 
assumptions made for our main Theorem 13. 1[ but we point out that, as can 
be seen from the proof given in Section for the validity of Lemma 14.21 it 
suffices to assume C^'^ regularity of dD and C°'^ regularity of the coefficients 
of the elasticity tensor C. 

We also need the following trace-type inequality. 

Lemma 4.6 (Trace-type inequality). Let Q be a bounded domain in M^, 
satisfying (12. 7p and (12. 9p . and let D be a simply connected open subset of 
fl satisfying (12. 8p and (12. lip . Let the fourth order tensor P be defined by 
fl2A8|) and satisfying (|2J9|) . (E^HD and ( K29^ . LetMe H-^dn,R'^) satisfy 
(I2.15p - (l2.16p ■ with S satisfying (I2.13p . Let w G H^{Q\D) be the unique weak 
solution of the problem (12.321) - (12. 37p . We have 

(4-7) l|M||^_.(^,^^,^<C||V^^IL.(,,^), 

where C > only depends on Mq, Mi, 6o and M. 

Remark 4.7. The proof of the above lemma can be obtained by adapting to 
the mixed problem the proof of the analogous Lemma 7.1 in [M-Ro-VeT] , 
which was established for the solutions to a Neumann problem. Indeed, the 
above result holds under weaker assumptions: it suffices that the coefficients 
of the elasticity tensor C are of class L°° and that dD is of class C^'^. 

Next lemma ensures regularity of the boundaries of a family of domains 
approximating the sets = {x E Q \ dist{x, dDi) < p}. 
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Lemma 4.8 (Regularized domains |Al-B-Ro- Ve] . Lemma 5.3). Let Q be a 

domain with boundary of class C^'^ with constants Pq, Mq and satisfying 
(12.71) . Let D be a simply connected open subset of Q with boundary of class 
C^'^ with constants po, Mq and satisfying (12. 8p . One can construct a family 
of regularized domains C Vt, for < h < apo, having boundary of class 
with constants po and Mq, such that 

(4.8) D C D''' C D''\ 0<hi< /i2, 

(4.9) 7o/i < dist(x, dD) < 71 /i, for every x e dD^, 

(4.10) \D'' \ D\ < j^Mipoh, 

(4.11) la^'^li <73Mipo, 
for every x G dD'^ there exists y G dD such that 

(4.12) \y - x\ = dist{x,dD), \i^{x) - u{y)\ < j^—, 

Po 

where i^i^x), z/(y) denote the outer unit normals to at x and to D at y 
respectively, and a, 7^, j = 0, 1, ...,4, and the ratios ^ and ^ only depend 
on Mq. Here | ■ |i denotes the 1- dimensional measure. 

Proof of Proposition \3.2[ The proof of this Proposition is rather technical 
and follows the lines of the proof of Proposition 4.2 in |M-Ro-Ve5] . which 
refers to the Neumann problem for the plate equation under the same gen- 
eral assumptions on the plate tensor and on the couple field. For this reason, 
in order to simplify the present exposition, we illustrate the main ideas un- 
derlying the proof of the first part of the result which is essentially the same 
as in the proof of Proposition 4.2 in |M-Ro-Ve5] . 

The main idea is that of applying the three spheres inequality ( 14. 2 p over 
a chain of disks centered at points of a path connecting x to any point 
y E {fl \ D)sp, thus estimating |V^typ in terms of ^-.^ |V^wp. In 

doing this, in order to reduce the number of iterations, which is responsible 
of the deterioration of the estimate, we exploit the regularity of the boundary 
components dVt and dD oiVL\D (in fact Lipschitz regularity suffices for this 
step), constructing suitable cones near the boundary, inside which one can 
apply the three spheres inequality over a chain of disks tangent to the cones. 

It is not restrictive to assume that po = 1- 
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By a covering argument, one obtains the following estimate, which is the 
analogous of (5.36) in |M-Ro- Ve5] : 
(4.13) 

-Ai-Bi log! 



B,(x) Jn\D \ L\d\^^'^ 



2„,,|2\ <5x 



vM' > I I — r — 1 yp<p 



where 6^^, < 6^ < 1, p, Bi and C are positive constants only depending on 
7, M, 6i and Mq, whereas Ai > only depends on 7, M, Si, Mq and Mi. 
Let us set 



By Holder inequality 
(4.15) 



l^'^lli2((n\I5)\(n\I5),.+,),) < |(^ \ ^) \ \ ^)(^+l)p|' l|V'^||'4((f,\i5)\(f,\i5) 



and by Sobolev inequality |Ad] 

we have 
(4.17) 

l|V^-lli^((.M5)v.\l5),,.,) < |(^^ \ I^) \ \ 11-11^ 

where C only depends on Mq, Mi. 

By interpolating the estimates fl2.38p and (14. 4p . we have 

(4.18) \\w\\^s^^^-^<C\\M\\m9n,R^), 

where C only depends on Mq, Mi, 7 and M. 
Moreover 

(4.19) |(^]\D)\(^]\D)(,+l),| <cp, 

with C only depending on Mq and Mi, see for details (A. 3) in |Al-Ro] . From 
dHZD, (Km\ and (Km> we have 



(4.20) / _ _ \V^w\^ <Cp^\\M\\l2 

J(Q\Z))\{Q\Z))(,+i)p 
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{9n,R2)) 



where C only depends on Mq, Mi, 7, M. 

By fli?Ti|) . and fHTj) we obtain that there exists p > 0, only de- 

pending on Mq, Ml, 60, J, M and F, such that for every p <p we have 



(4.21) 



/(n\D)(,+i)^ 



1 

> -. 

- 2 



By (KT^ . flOT]) and (iT]), we have that for p < min/ p, p}, 



(4.22) 



V^w\^ > (Cp 



Bp{x) 



iMf 1 



where C > only depends on 7, M, 61, Mq, Mi and 5o- Let us take p < C. 

Noticing that | logp| < ^, for every p > 0, and that p < 1, by straightforward 

computations we obtain that (13.51) holds with A = 3exp(j4i| log^^^l), B = 

I \ogS^\Bi + 1 for every p < p* with p* = mm{p, C}, p* only depending 
on 7, M, 61, Mo, Ml, 6q, and F. _ 

If p > p* and X e {Q\ D)sp C (^7 \ -D)sp., then, by applying (13. 5 p with 
p = p*, we have 



(4.23) 



B„,(x) 



\V^w\^>C*\\Mf 1 

' I — II 



where C* only depends on 7, M, 5i, Mq, Mi, and F. 
Since x G (f2 \ -D)sp, we have that 

(4.24) diam(fi > 2sp, 
and, on the other hand, 

(4.25) diam(0\D) < C2, 
with C2 only depending on Mq and Mi, so that 

(4.26) 

By fl4:23|) and (1^26]) . we have 

(4.27) / \V^w\^ > 

Jb 



2s 1 

C'2-~p 



C 



'Bp(x) 



exp 



A 



B 



Mf 1 



with C = C* exp 



□ 
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Proof of Proposition \3.3[ Let G be the connected component of Q\ [Di U D2) 
such that S C dG. Let us prove dSSD-M, the proof of ([X7D-(IS3D being 
analogous. 

Let 'g ^ A the affine function reahzing the minimum in (13. ip and let us 

set 

(4.28) w = wi — W2 — g. 

Step 1. Let ^ = 0. 

In this case w = Wi — W2- 

It is not restrictive to assume e < Po\\M\\ 1 Jl, where Jl, < Jl < 

is a constant only depending on 7, M, 5i, Mq, Mi, which will be chosen 
later on. 

Let ■& = min < a, 7^, ^" >, where /iq, < /iq < 1, only depending 

[ ''Ti 47o-\/l+M,f J 

on Mo, is such that Qh is connected for every h < ho (see Prop. 5.5 in 
|Al-R-Ro-Ve] ). and where a, 70, 71 have been introduced in Lemma H78l Let 
p = -dpo and let p <p- 

Let us denote by Vp the connected component of i7 \ (D^ U which 
contains dQ, where are the regularized domains introduced in Lemma 
01 We have 

(4.29) D2\d;c ({d{ \d[)\g)u {{n \ %) \ , 

(4.30) a((fi\v;,)\D^) =f?uf^, 

where = dD!^ n dVp and f? C 9/)^ 
We have 

(4.31) / _ivvr</ ivvr+ / _ivvr. 

By (iSl) and flirTOj) we have 

(4.32) / \V'w,\'<Cpl\\M\\l_. 

J{D^\Dl)\G " 2(an,R2)po 

with C > only depending on 7, M, Mq and Mi. 



20 



By applying the divergence theorem and by fl2.29p we have 



where 
(4.34) 



(4.33) / _|VVr<7"W B2{wi)wi,, + B2{wi)wi,s + B:i{wi)wi- 



B2{wi) = {FV^wi)u ■ u, 



(4.35) 



(4.36) 



Bsiwi) = -div (FV^wi) ■ ly, 



with 1/ denoting the unit outer normal to {Q\Vp)\Di, and s an arclength with 
associated parametrization (p{s), such that (p'{s) = T{ip{s)), where r = e^xu. 

Let X e tl. By (gj]), dist(x,aDi) < 7ip. Let y e dDi such that 
\y — x\ = dist(2:,9Di) < 7ip. Since 7ip < |po and by fl2.8p we have that 
a; G (fi \ D7) \ Uf°/^ and (USD applies. Recalling also that Wi = 0, Vwi = 
on dDi, we have that 



(4.37) 



\w,ix)\ = \w^{x)-w,{y)\<Cpl\\M\\^^ 



3 (9n,M2) 



(4.38) Mx)\ = \{Vw,{x) - VwM) ■ < Cpo\\M\\^_^^^^^^^^-^^ 



(4.39) \w,,s{x)\ = \{Vw^{x) - Vw^iy)) ■ t{x)\ < Cpo\\M\\^_^^^^^^^^-^^ 

where C only depends on 7, M, Mq and Mi. 
From (glD, dOHD, fl4:39|) and f HTT|) we have 



(4.40) 



B2{Wi)Wi^y + B2{Wi)Wi^, 



where C only depends on 7, M, Mq and Mi. 
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(4.41) 



B'i{wi)wi 



< Cpo||53(Wl)||i2(fP)||w^l||i2(fP) < 



< ^ll''^i|lH3(f5')ll'"^i|lL°°(ff) ^ ^ll'"^ill/f4((n\Dr)\w«)/«)ll^illi°°{f5') - 



Po 



where C only depends on 7, M, Mq and Mi. 

Let X e f^. By (|M]), dist(x,aD2) < 7ip. Let y e ^ch that 
ly— xl = dist(x, 002) < 7ip. Since 7ip < |po we have that x E {VI\D2)\Up°^^ 
and fl4.6p apphes. RecaUing also that W2 = 0, Vw2 = on dD2, we have that 

(4.42) \w^{x)\ = \w{x)+W2{x)-W2{y)\ < + CPoll^ll^-i(9^,j,2) A 



(4.43) 

\wi {x)+{Vw2{x) -Vw2{y))-iy{x)\ < \Vw{x)\+Cpo\\M\\^^i 



P 



(4.44) 

\wi,s{x)\ = \w,s{x)+{^W2{x) - Vw2{y))-T{x)\ < |V«;(x)|+Cpo||M|| 

where C only depends on 7, M, Mq and Mi. 
From (gSl) fOHjl . flOill and fl4TT]l we have 



(4.45) 



B2{Wi)Wi^y + S2(tyi)tfl,, 



<Cpl\\Mf 1 — +Cpo||M|| 1 ,max|\ 



where C only depends on 7, M, Mq and Mi. 

Moreover, by fl4.42p and (14. 5p . we have, similarly to above. 



(4.46) 



^ 2 



53(u7l)«^l 



<Cp^||Mf 1 — + C||M|| 1 ,max|zi;| 



where C only depends on 7, M, Mq and Mi. 

By fiOTjl -f iOajl . flCTjl . fl44T]l . fH:i5|) and flCTj) we have 
(4.47) 
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where lkllci(eyp\on) = '^^^avp\dn \M + Po'^^^dVP\dn and C only de- 

pends on 7, M, Mo and Mi. 

Let us notice that w G H'^iG) satisfies 



(4.48) 
(4.49) 
(4.50) 
(4.51) 



f div(div(PV2^i;)) = 0, in G, 

{W'^w)n ■ = 0, on S, 

div(PV2«;) ■ n + ((PV2«;)n ■ r)„ = 0, on S, 



By the regularity of w near E (see Theorem 15.11) and by standard interpo- 
lation inequalities, we can apply to w the stability estimate for the Cauchy 
problem contained in |M-Ro-Ve4l Theorem 3.8]. Therefore, noticing that, by 

(I238D, \\w\\l2(g) < \\wi\\LHn\-Dl) + \\w2\\L^(^n\-m) < ^PoW^Wn-^idnM^y "^^^^ 
only depending on 7, Mq, Mi, we have 

(4-52) ikii..,«^,^,P.,n„, < c.v?'-''iiM|i;;:,„„,.,,, 

where C > 0, r G (0, 1) only depend on 7, M, Mq and Mi. 
Let 

(4.53) ^0 = - jn, 

where n denotes as usual the outer unit normal to dfl at Pq, and notice that 

PQ PQ 

By (1132D, 



(4.54) \\w\\ I X <Ce>f-^)||Mf-\ 

L2 B pn (^0) 

Let us consider the set fl fi_poJHj_- By the choice of p and recalling that 

/lo < 1, we have by straightforward computations that this set is connected 
and contains zq. If a; G y fl ^2 ppfap then, by the choice of p, dist(a;, dQ) > 

> 7oP and by dist(x, dDi U ^Da) > 7op. Therefore any disk of 
center in fl ^2 pphp and radius 70P is contained in G. 

Let X be any point in dVp\dfl and let 7 be a path in V^nQ pphp joining 
x to 2;o. Let us define {xi}, i = l,...,s, as follows: xi = zq, Xi+i = ^{ti), 
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where = max{t s. t. \j(t) — Xi\ = if \xi — x| > otherwise let 

i = s and stop the process. By construction, the disks BjoR{xi) are pairwise 
disjoint, — Xi\ = for i = I, s — 1, l^^ — x\ < Hence we have 

s< S (^^y, with S = only depending on Mq and Mi. 

An iterated application of the three spheres inequality (14. 2 p to w over the 
disks of center and radii ri = r2 = = 7op, gives that for every 

p, < p < p, 

/ \ 5" 

2 \ I / l„.,|2 



(4.55) / \w\' <C[ \w\' / \w 

JBji^ix) \Jg / WB7^(20) 

where S, < 6 < 1, C > 1, only depend on 7, M and 5i. 

Since Bjop{zo) C B po (-^n). by l^iM), fli35|) and flOHjl we have 

(4.56) / |«;r<Cp^i|Mf e^^^\ 

where 
(4.57) 



Pnll^ll 1. 



Let us recall the following interpolation inequality 



,4 1 1 , , 

2 \ Irr 12 I / / I |2 



(4.58) ||t;|U^(B.) < I ( / ) |V^;|1.(5^) + - 

which holds for any function v G W^'°°{Bt). 

By applying (14.581) to w in i?70p(2;) and by (I4.56P and (14.60 . we have 

By the following interpolation inequality (see for instance |G-T] ) 

1 1 

(4-60) IMcHB,) < (^11^11^0(50 



we have, recalling (14. 6p . that 



(4.61) IkllcHi^^M) < ^Pol|M||^..(,,^,,) ] 
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Po\^ -111 



so that 

(4.62) \Hy,oy.,,,, < Cpl\\M\\,-i,,,^,., (^) ^ 
From fl4.47p and fl4.62p we have that, for every p < p, 

(4.63) / I < cpiwmu^,^^,, + ^e^) . 

Let us set p = exp | — ^exp ^ ^-^I^Qg*^! j |^ ^ = niin {/i,exp (— ^)}. We have 
that p, p < e~^, only depends on 7, M, 5i, Mq and Mi. Let e < p, and let 



(4.64) p(e~) = po 



2g|log(5| 

log I log 6 4 



Since p(e) is increasing in (0, e ^) and since p(/i) < p{p) = po9 = p, we can 
apply inequality ( I4.63P with p = p(e), obtaining 

(4.65) /^^^_|VV,P<Cp?||M||^_j^^^,,^(log|log?i|)-', 

where C only depends on 7, M, Mq and Mi. 

Since e < exp(— ^f), we have that log j > —\ log | loge|, so that 

(4.66) log |log £4 1 > - log |log e| . 
From fl4.65p and fl4.66p we have 

(4.67) X^^_|VV,r^<P?l|M||^_.^^^,_^.(E), 
with 

(4.68) uj{t) = C(log|logt|)"^ for every t < 

where C > is a constant only depending on 7, M, Si, Mq and Mi. 
Step 2. Let g ^ 0. 
Let 

(4.69) P = p{~e), 

where p{e) is given by fl4.64p . 
We have two cases: 



25 



I) aDfnr^ = 0; 

II) aDfnf§^0. 

When I) holds, there are three possible subcases: 
la) b{r}D^ = 0, 

lb) b{ c , 

Ic) Dl c bl 

In case la) we have that {n\Vp)\b{ = and, therefore, d {{VL \ Vp) \ b(^ = 

dbl whereas in case lb) (fi \ Vp) \b{ = D^\b{. 

For both cases, by applying the divergence theorem to wi in and in 
D2 \ Di respectively and by taking into account f l2.37p . we have 
(4.70) 

/ (div(PVV) ■ ((PVV)i^ ■ t)„) £/-((PVV)z^-z^)£/,i. = 0, for every G 

with V denoting the unit outer normal to and s an arclength with asso- 
ciated parametrization v'(s), such that </?'(s) = r((/j(s)), where r = 63 x z/. 

By applying the estimates of continuation from Cauchy data (14.621) ob- 
tained in the above step to the function w defined by (I4.28p . we have 

1 

By recalling that = 0, Vwj = on dDi, i = 1,2, and by arguing similarly 
to Step 1, we have, for both cases, that 

(4.72) / _|W|'< / _|VV|'< 

Jd2\Di Jd^\Di 

< / (div(PV^zi;i) ■ + {{FV^Wi)u ■t),s)wi- {{FV^Wi)u ■ i^)wi ^ = 
JdDl 

= 7"^ [ (diviFV^Wi) ■ V + ((PV2u;i)z^ ■r),,)w2- ((PV^w;i)z^ ■ v)w2 u+ 

+ 7~W _ (div(PV^u;i) ■ v + ((PV^ )w - {{FV^Wi)iy ■ iy)w,^< 

JdDl 

<Cpl\\Mf 1 ( I + P^i¥'] = plWMf , uCe). 

- '^Oll "H-3(9f7,R2) p J "H-^{dn,R") 

In case Ic), by using (14.101) . we have 

(4.73) |I^2\^| < \bl\D,\<Cpop, 
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with C only depending on Mq, Mi. By fl4.73p and by (14.61) . we have 
(4.74) /^^^_|VV,|^<Cp„p||M||^_.,^„„,,, 

with C only depending on 7, M, Mq, Mi, so that the thesis follows. 
Let us consider now case II). In view of the above arguments, 

(4-75) M\c^iav,\an) < ^Pol|M||^_i^^^^^,^a;(e~), 

where C > only depends on 7, M, 5i, Mq, Mi. 
Let z e 9Z)f n f We have that 

(4.76) \w{z)\ + po|V«;(z)| < Cp^||M||^_, ^^^^^^^c^(e), 

where C > only depends on 7, M, ^i, Mq, Mi. On the other hand, by 
using the homogeneous Dirichlet conditions on the boundaries of the rigid 
inclusions, 
(4.77) 

|^.(z)| +Po|V^.(.)| < Cpl\\M\\,A,enmi ^ ^^ol|M||^.. (,^^^,^a;(e), 
where C > only depends on 7, M, 5i, Mq, Mi. By fHT5]) - fHT7j) . we get 
(4.78) 1^(^)1 + po|V^(^)| < Cpl\\M\\^,r^^^^^^^uii). 

Let ^(x) = axi + bx2 + c. Then, by (14. 78 p . we have 
(4-79) |a|<Cpo||M||^-i(,^,,.)C.(e~), 

(4.80) \b\ < Cpo\\M\\^-i,,^^^. fi^, 

(4.81) |c| < 1^(^)1 + \a\\z\ + \b\\z\ < Cpl\\M\\^_i^^^^^^^u{~e), 

where we are assuming for simplicity that the origin belongs to Q. It follows 
that 

(4.82) mcH^<Cpl\m^-^,,,^^^.^<e)- 

By repeating the arguments of Step 1 ioi w = wi — W2 — g, we have the addi- 
tional term |^| which is controlled by fl4.82p . and the thesis follows similarly. 

The proof of the better rate of convergence (13.90 . under Lipschitz regular- 
ity condition on the connected components of dG can be obtained by merging 
the geometrical construction illustrated in |M-Rolj and the arguments seen 
above. 

□ 
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5 Proof of Lemma 



12 



and Lemma 



as 



Let us denote by 

(5.1) B+ = {(2/1, y2) G M'l + yl < a\ > 0} 
the hemidisk of radius a, ex > 0, and let 

(5.2) r, = {(7/1,2/2) GM^I -a<y^<a, y2 = 0}, 

(5.3) Tt = dB^\T,. 
Moreover, let 



(5.4) H^AB^) = ^ge H\B+)\ ^ = 0, ^ = on 

Proof of Lemma \4.^ It is not restrictive to assume po = 1. By the regularity 
of dQ and dD, and by (12. 8p . we can construct a finite collection of open sets 
fio, {^j}-Li, such that 

(5.5) {vij r]n)c{n\D), j = i, a^, {n- r]n)c{n\D), % = i, a^', 

(5.6) fijnfi; = 0, i = l,...,A^', i = l,...,iV, 

(5.7) 1] \ D = fio u (uf=i7^^^(i?t)) u (ufJir(;J(i?t 



(5.8) r]oc(r]\D) 



(5o 1 



where 5o only depends on Mq. Here, 7(j), j = 1, A^, is a homeomorphism of 
C^'^ class which maps QjCiQ into Bf, QjddQ into Fi and dQj n{Q\D) into 
r^. Similarly, T\i) is an homeomorphism of C^'^ class which maps fl f2 
into B^, TVi n into Ti and dn[ n (il \ D) into T^, i = 1, A^'. It can be 
shown that every mapping T(j), T\i), i = j, N , i = 1, N', can be chosen 
such that the Jacobian of the transformation is identically equal to one, see 
|Ag| (p. 129). By the regularity of dQ and dD and by (12. 7p . the numbers 
A^, N' are controlled by a constant only depending on Mq and Mi. 

By covering Qq with a finite number of spheres contained in Q\ D and 
using local interior regularity results (see, for instance, jM-Ro- VeT] . Theorem 
8.3), we have that w G H^{Qo) and 

(5-9) 1^11^3(^0) < C\\w\\H^(n\D), 
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where the constant C > only depends on Mi, ||P||c70,i(n) 7. 

Let us fix J, 1 < j < A^, and let us show that an analogous estimate holds 
true for Qj fl Q near the boundary dQ where the Neumann data M is given. 

The function w, solution to fl2.32p - fl2.37p . satisfies 



(5.10) 



We define 
(5.11) 



M„(/j,„ + {Mr),s^) ds, 

for every ip G //Jn,nn(^- ^ 



y = T(j){x), yeBf, 



(5.12) 



(5.13) 



u{y) = w{T^j^\y)). 



Then, changing the variables in fl5.10p . the function u belongs to H'^{B^ 
and satisfies 



(5.14) 
where 



a+{u,^p) = l+ii/j), for every i/j G H^+{B^), 



(5.15) a+(-u, Ip) = 



det 



577 



(i) 



dx 



-1 



dy, 



(5.16) l+{^) = - 



p, \ OX \ ox 



dT, 



dx 



n 



+ 



0'), 



dx 



dT(. 



dx 



d^, 



with 
(5.17) 



Sk 



dXr 



R. 



■ksr 
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(5.18) Qlkmn ^ ^ -Pijrs^lsSkrSmiSnj y 

i,j,r,s=l 

2 

(5.19) -^Ikn ^ ^ PijrsSlsSli^rRnij y 

i,j,r,s=l 

2 

(5.20) Bj^jyiji ^ ^ PijrsSmiSnj-Rksr y 

id,r,s=l 

2 

(5.21) D]f^ ^ ^ PijrsRksrRnij 

i,j,r,s=l 

and 

(5.22) ^(y) = ^(7^^i(y)), V^e//2^(i?+), 

(5.23) >l„(y) = M„(7^^i(y)), Mriy) = M;(7^^'(y)) 

and where p is the unit outer normal to B^. By fl5.18p and fl2.19p . 02.291) . 
the components of Q satisfy the symmetry conditions 

(5.24) Qafi-yS = Q^Safi = Q-ySfia, «, = 1) 2, in , 

and the strong convexity condition 

(5.25) QA-A>7o|v4|^ in5+, 

for every 2x2 symmetric matrix A, where 70, 70 > 0, is a constant only 
depending on 7 and Mq. 

By the regularity assumptions on P, the tensors Q, A, B, D belong to 
C°'i(S+). Moreover, G H^T^) and Mr,^ G H-^T^). 

Now, we use the following regularity result up to the boundary for the 
solution u to the problem (15.141) . which is proved in the Appendix. 

Theorem 5.1 (Boundary regularity for non- homogeneous Neumann condi- 
tions). Let u G H'^{Bi) defined by (15.131) be the solution to (15.141) . where 
the tensors Q, A, M, D and the couple field Ai are defined as above. Then, 
u E ("^i ) have 

(5-26) \\U\\ / ^\ < C (\\M\\l^^ + \\u\\rr2lF!+) 

V J II II ^3 1^5+ j - V" "H^iFi) II ll^^(Si) 

where C > only depends on Mq, ||P||c'0,i(q) and 7. 
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By applying the homeomorphism 7(j) to (I5.26p we have 

(5.27) / / \\<c(\\M\\ i, + Ikllffamnm 
^ ^ II "m(T^:l(^B+Jj- V" iiH2(n,n9n) ^ II nn (u,nu) 

where the constant C > only depends on Mq, 7 and ||P||(70,i(q)- 

We now derive an estimate analogous to fl5.27p near the boundary of the 
rigid inclusion D. Let us fix z, 1 < z < A^'. The function w, solution to 
(^M-<^M, satisfies 

(5.28) / PV^w ■ VV = 0, for every (f e H^{Q'i n Q). 
Jn'-nn 

Then, by introducing the transformation 

(5.29) y = r\,){x), yeBt, 

(5.30) x = T'i^{y), xen'^nn, 



(5.31) viy) = w{ri^{y)), 

and changing the variables in fl5.28p . the function v belongs to (Bi) and 
satisfies 

(5.32) a'^{v, i)) = 0, for every V e H^iBf), 
where 



(5.33) a'^iv,!!;) 



Jb+ 



dx 



-1 



and the tensors Q', A', B', D' and the function -0 are defined as in fl5.17p - 
fl5.2ip and (15.220 . respectively, with T(j) replaced by T'(j). We note that the 
tensor Q' satisfies the conditions fl5.24p . fl5.25l) and all the tensors Q', A', B', 
D' belong to C^'\Bf). 

To this point we make use of the following regularity result up to the 
boundary dD for the solution v to the problem fl5.32p . 
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Theorem 5.2 (Boundary regularity for homogeneous Dirichlet conditions). 

Let V G H^_^{B^) defined by (15.311) be the solution to (15.321) . where the tensors 

Q' , A' , M' , D' are defined as above. Then, v G ("^^ ) have 
(5.34) ll'^llHafBl) ^'^MnHBty 



where C > only depends on Mq, ||P||(70,i(T7) o,nd 7. 

The proof of Theorem 15.21 foUows the same guidehnes of the proof of 
Theorem 15.11 see also |Ag] for details. 

By applying the homeomorphism T'(i) to (15.341) we have 

(5-35) 1^11^3/ <C\\w\\H^n'^nn), 



where the constant C > only depends on Mq, 7 and \\ir\\QO,i(Qy 

Finally, estimate (jOl) follows from ([52]), (^M>, flHTTTl) . §3B and from 
the estimate □ 

The proof of Lemma follows from the following local version of the H"^- 
regularity near a boundary with homogeneous Dirichlet data for the solution 
to the problem (I5.32p . 

Theorem 5.3. Let v G H^^{Bl) defined by fICT]) be the solution to ([532]). 

Then, v G H^{Bt) and we have 
2 



where C > only depends on Mq, ||P||(7i,i(n) ^'^^ 7- 

Proof. By Theorem 15.2] the function v G H^{Bt) satisfies the estimate 

2 

(I5.34p . e.g. \\v\\ ^/ \ < C\\v\\^2i^B+Y ^'^^^^ C > only depends on M 



co.i(n) and 7, 



0, 



In the first step of the proof we estimate the tangential derivative 
By the weak formulation of the problem (I5.32p it can be shown that the 
function v^a = ^ £ l"^^)' ~ satisfies the following equation 

(5.37) a'^v,^, ij) = l'^4^), for every i; G H^^ (st) , 
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where 



(5.38) a'^{v^^, i^)= I (Q'V^^;,^ ■ V^V' + A'V^t;,^ ■ 

J 



det 



5x 



-1 



dy 



(5.39) /UW= / (div (Q;„ ■ V'^ - A',«V'^ ■ 



det 



dx 



dy. 



By the regularity assumptions on P and by (15.341) we have 
(5.40) \lUm<C\\v\\ 



a = 1, 2, where the constant C > only depends on Mq, and 7. 

To this point, the arguments used in the proof of (I6.35P for the Neumann 
case can be adapted to the Dirichlet boundary condition case to obtain 



(5.41) 



oyi 



L^lB+ 



<C\\v\\ 



H'2(B+)^ 



a = 1,2, where the constant C > only depends on Mq, ||P||ci,i(n) aiid 7. 



The estimate of the normal derivative f a = 1,2, can be obtained 



by adapting the proof of (I6.44p . We obtain 

„ d 



(5.42) 



dy2 



V'v. 



'L2lB+ 



<C\\v\ 



H'^{B+)^ 



a = 1,2, where the contant C > only depends on Mq, ||Pilc'i,i(Q) and 7. By 
fl5:iT]) and ([532D we have the thesis. □ 



6 Appendix 



In this Appendix we prove Theorem 15.11 
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Proof of Theorem \5.1[ As a first step of tlie proof, let us estimate tlie partial 
derivative -^V'^u in tlie direction parallel to the fiat boundary Fi of B^. 

Let I? e d^(]R^) be a function such that < 'd{y) < 1 for every y G M^, 
7? = 1 in 5p, t9 = in \ and < C, k = 1, ...,4, where p = |, 

^0 = |(p + 1) = I and C is an absolute constant. 

For every function ip G H^+{B^), we still denote by -0 G if^(]Ri) its 

1 

extension to obtained by taking -0 = in \ Sj^. 

Let s be a real number different from zero. The difference operator in the 
yi-direction is defined by 

(6.1) (r,„/)(,) = f(y^^^')-ny\ 

s 

for any function /. In what follows we shall assume that |s| < jg- Let us 

notice that if w G H^iBl), then n^si^u) G Hl+{Bt). 

^ 1 

We evaluate the bilinear form a+(-,-?/') with u replaced by T\^a{'du) and 

ij G Hl+{Bt). Since 
^ 1 

(6.2) V\t^,s{^u)) = ri,,(V'=(^«)), in S+, A; = 1,2, 
we have 



(6.3) a+(Ti,,(i9M),V^) 



/ (Qri,,(V'(^9M)) ■ + Ati,,(V2(^m)) ■ VV^+ 
+lri,,(V(^?M)) ■ V=^^ + Dti,,(V(^9m)) ■ Vt/-) (iy. 



where we have taken into account that det = 1. 

Let us consider the leading term of (16.31) . We elaborate the corresponding 
expression by moving the difference operator from du to dif). We have 



(6.4) / Qri,,(V'(i9n)) ■ VV = / QTi,,(W'ix) ■ V'^+ 

lB+ Jb+ 



Jb+ Jb+ 
The last two terms on the right hand side of (16. 4p can be estimated as follows 

(6-5) I/2I <C\\u\\H2^^B+)\\^'^'^\\L^{B+)^ 

(6.6) I/3I < C'||n||^i(^+)||V^V'llL2(B+), 



34 



where the constant C > only depends on ||IP||ioo(^+) and Mq. 
The term Ii can be written as 



(6.7) h = [ 



Jb+ ^ 



where 
(6.8) 



\n\<C\\V'u\\LHBt)\\^'nL^^BtV 

with C > constant only depending on 

The remaining term I[ can be elaborated as follows 



(6.9) l[= I n,. 



'm) ■ (ri,_,V'V^) 



Q(V\)-V2(i?ri,_,V^) + 



+ 



'u 



The last two terms of fl6.9p can be estimated as follows: 



(6.10) 



■n ■ ((V2t9)ri,_,^) 



<C\\WML^i,Bt)\\^nLHBtY 



(6.11) 



< 



< C||V'm| 



L^B+)\ 



\LHB+)^ 



where the constant C > only depends on ||P||ioo(^+) and Mq. 

Therefore, by ([621), (ESD, dHSD, ([63]), (16^01) and fl6llD . the left hand 
side of (16.41) can be written as 



(6.12) [ Qri,,(V2(i9M))-VV = - [ 
Jb+ Jb+ 



where, by the Poincare inequality on H^+{Bi), we have 



(6.13) 



< C'll^i|l/f2(B+)l|VVllL2(iJ+), 
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where the constant C > only depends on ||P||^(, -^^-^^ and Mq. 

We can estimate the remaining terms appearing on the right hand side 
of (16. 3p similarly. Concerning the term involving the tensor A, for instance, 
we have 

(6.14) / Ari,(V^(i9M)) • = / KTi,{dV'^u) 
Jb+ ' Jb+ 

+ Ari,(V^9®VM + VM®V^9)-VV^+ / An ,(mV¥)-V^ = Ji + Ja + ^s, 

Jb+ ^ Jb+ 

where, by the Poincare inequality on H^+{Bf) we have 

^ 1 

(6.15) I Ja + J3I < C\\u\\jj2(^B+)\\^i^\\LHB+)^ 

where the constant C > only depends on ||P||ioo(5+) and Mq. The term Ji 
can be elaborated as it was done for Ii, obtaining 



(6.16) Ji= TisiA{^V'u))-V^- 
iBt ' 



- [ {n,sA){^V\){y + sei) ■ = J( + Jf, 
Jb+ 



where 
(6.17) 



\Jn<C\\VMLHBt)\\^i^\\LHBt) 



with C > constant only depending on ||P|Lo ^o- Concerning the 

term J[ we have 

(6.18) J[= [ Ti ,(A(?9V\)) • = - / AV^u ■ (i^n _,V^) = 
Jb+ ' Jb+ 



A{V^u)-V{^n.,ij)+ / AV^M- V79(ri_,V), 
Bt ' JBt 



where the last term of (16.181) can be estimated as 



(6.19) 



Jb+ 



<c\\ liv^ll 



with C > constant only depending on ||P||ioo(B+) and Mq. 
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Therefore, by ( KI^ . flKTTj) . fETTSj) and flCTjl . the left hand side of IKT^ 
can be written as 



(6.20) 



Jb+ ' Jb+ 



where 

(6.21) I 
with a constant C > only depending on 



< C'lkllH2(B+)IIWllL2{B+). 



By similar procedure we obtain 



co,i(^) and Mq. 



(6.22) 


[ Bri,,(V(i?M)) ■ 


= - 


1 BVu-V^{dTi.si') + 




Jb+ 




Jb+ 


(6.23) 


[ Dri,,(V(^9n)) 


■ = - 


/ DVm ■ V(t9ri -.V) + ' 










where 








(6.24) 


|rB| < C\ 


'^\\h^{b+) 




(6.25) 


\ro\ < C 


\W\\h'^{b+) 


I|vV'IIl2(b+), 



where the constant C > only depends on ||P||^o,i(;^) Mq. Collecting 
the above results and by using the Poincare inequality on H^+{B^), by (16.31) . 
IKT^ . (Km . we have 

(6.26) a+^Ti^si^u),^;) = -a+{u,{}Ti_si') + r+, 

where 



(6.27) 



< C'||m||^2(b+)||VV|Il2(b+ 



with C > constant only depending on ||P|| ^ -— and Mq. Since tp E 
H^+{B^), the function dTi^-si^ G H\{B^) is a test function and then, by 
the weak formulation of the problem fl5.14p . we have 

(6.28) a+(M,^9Ti,_,V^) = /+(i9Ti,_,V^) = 



M 



9T{j) f dTjj) 
^ dx \ dx 



+Mr. 



9%) 



dTi 



(i) 



dx 



dx 



n 



+ 



U). 



dx 



-1 
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By trace inequalities and Poincare inequality we have 
(6.29) |/+(i9ri,_,^)| < 



'//2(ri); 



where the constant C > only depends on Mq. By fl6.26p - (l6.29p we have 

72„/,|| 

■)' 



(6.30) a+(Ti,,(^M),^) < Ci (\\M\\^. + ||V^^||^,(^^ 



for every ip G H^+{Bi), where the constant Ci > only depends on \\co,i(b+) 
and Mq. 

Let ip ^ and let us estimate from below a+i^ip,!/)). For every 

^ 1 

e > and for every tjj G H^+{B^) we have 

i 1 



(6.31) 



< 



1 



<C(e||VVl|i.(^+)+(l + -)||VV'||J,(^.) 



co,i(^) ^0- Therefore, 



where the constant C > only depends on 
by the strong convexity of Q and choosing e small enough in fl6.3ip we have 

72 -,||2 IIV7„;.I|2 



(6.32) a+{^P,^) > C^H V^H^.^^^) - C3II V^H^.^^^ 



)' 



where C2 > 0, C3 > are constants only depending on ||P||^g 7 and 
Mq. Now, by taking t/j = Ti^si^u) in (16.300 and (I6.32p . we obtain 

(6.33) C2\\V\n,,{du))\\l,^^^^ < C,\\V{T,,s{m\\l.^^^^^+ 

+ C, + 11^11^2(5+)) ||V2(ri,,(^^z))||^.(5+), 



where the constants Cj > 0, z = 1,2,3, only depend on ||P||^g 7 and 

Mq. Therefore, recalling that || V(ri^s('(?'ii))||L2(_B+) — c||m||j:^2(5+), where c is 
an absolute constant, and by Poincare inequality we have 



(6.34) 



|V'(ri,,(^9n))| 



<C[\\M\\i,^ . + \\u 
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where the constant C > only depends on ||P||^o 7 and Mq. Taking 
the hmit as s — )■ we finally have 

d_ 

^2/1 L2(B+) 



(6.35) 



where the constant C > only depends on ||P||^(, 7 and Mq. 

To obtain an analogous estimate for the normal derivative ^V^m, we 
shall use the following Lemma. 

Lemma 6.1. ( lAgl Lemma 9.3/^ Assume that the function g G L'^{B^) has 
weak tangential derivative ^ G L'^{B^) and that there exist a constant K, 
K > such that 



(6.36) 



dy'2 



< mWmiB^), for every ^ G C^{Bt). 



Then, for every p < a, g E H^{B^) and 
(6.37) ll^ll <C 1^+11^11x2(5+)+^ 



dg 




dyi 


LHB+)) 



where C > only depends on ^. 

Throughout this part let -0 G C^{B^^). The arguments above show that, 
without loss of generality, we may assume that A = B = 13 
we can write 



(6.38) a+{u,i:) 



= B = D = 0. Therefore, 

{Q22'y5U,^5 "0,22 +2(5l27<5li,75 'V^,12 +Qll-i5U,^5 '0,ii ) 



and, by integration by parts, we have 
(6.39) 

Q227<5M,7<5 = a+{u, i>) + 2 {Ql2^5,l U,-yS tp,2 +Ql2-,5U,^51 4',2 ) + 



+ / iQiij5aU,.ysil'a+Qn-ysu,^sii'a) ■ 
Since ifj G C^{B'^^), we have aj^{u,ip) = 0. Let us define 

2 

(6.40) 9 = ^ Q22-,5U,^5 ■ 

7,5=1 
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Clearly g e L\B+) and ^ e L'^{B+) by (K^ . By flCTj) and by estimate 
fl6.35p we have 



with C > constant only depending on ||P||^o 7 and Mq, so that, in 

5+ , the function g satisfies the conditions of Lemma [6TT1 Hence, g G H^{B^) 
and we have 

(6.42) Mhhb^) < C + , 

where the constant C > only depends on ||P||^j, j^-^^, 7 and Mq. By the 
ellipticity of Q, by ( K^ and it follows that 

(6.43) U,22 = Q2"222 ^ - Yl Q22^SU,^S G ^^'(5;), 

\ (7,5)^(2,2) / 

and therefore 

(6.44) \\U,222 \\mB+) < C (llAlll^l^^^^ + , 

where the constant C > only depends on ||P||po,i(B+)' ^ 
By flOHjl and flCTjl we have 

(6.45) < C (llA^ll^i^^^^ + , 

where the constant C > only depends on ||P||^q -^j,-^^, 7 and Mq. This 
completes the proof. □ 



(6.41) 
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